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ABSTRACT 

We study the K-correction for the case of emission hnes formed in the X-ray iUuminated atmosphere 
of a Roche lobe filling star. We compute the K-correction as function of the mass ratio q and the disc 
flaring angle a using a compact binary code where the companion's Roche lobe is divided into 10^ 
resolution elements. We also study the effect of the inclination angle in the results. We apply our model 
to the case of the neutron star low-mass X-ray binary X1822-371 (V691 CrA), where a K-emission 
velocity Kem = 300±8 km s~^ has been measured by Casares et al. (2003). Our numerical results, 
combined with previous determination of system parameters, yields I.6IM0 < Mjvs < 2.32M© and 
OAAMq < M2 < O.56M0 for the two binary components(i. e. 0.24 <q< 0.27), which provide a 
compelling evidence for a massive neutron star in this system. We also discuss the implications of 
these masses into the evolutionary history of the binary. 

Subject headings: accretion, accretion disks - binaries: close - stars: individual: X1822-371 - X- 
rays:stars 
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1. INTRODUCTION 

Low mass X-ray binaries (LMXBs) are interacting 
binaries where a low mass star transfers matter onto a 
neutron star (NS) or black hole (e.g. Charles & Coe 



2004). There are ~200 bright (Lx 
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LMXBs in the Galaxy and most of them harbour 
neutron stars as implied by the detection of X-ray 
bursts/pulsations. About ~20 of these have period 
determinations but only in very few cases (e.g. Cyg X-2: 
Casares et al. 1998, Orosz & Kuulkers 1999 and V395 
Car: Shahbaz et al. 2004, Jonker et al. 2005), a full 
orbital solution and mass determination exists thanks to 
the spectroscopic detection of the giant companion star. 
For the others, the companion is completely swamped 
by reprocessed light from the accretion flow and this 
has hampered detailed knowledge of their dynamical 
properties. 

X1822-371 (V691 CrA) is one of the brightest LMXBs 
in the optical, and also the prototypical accretion disc 
corona (ADC) source due to its low L^/Lopt ~20 ratio 
(it is usually ~500-1000). In ADC sources, the combi- 
nation of high inclination and a thick disc obscures the 
central source, and only X-rays scattered from material 
above and below the disc can reach the observer (White 
& Holt 1982). Thanks to eclipses of the accretion disc 
by the companion star on a 5.57h orbital period (Hellier 
& Mason 1989; hereafter HM89) an accurate determina- 
tion of the inclination angle exists, i=82.5 ±1.5 (Heinz & 
Nowak 2001; hereafter HNOl). The source also displays 
0.59s X-ray pulses (Jonker & Van der Klis 2001) which 
enables an extremely precise determination of the neu- 
tron star radial velocity semiamplitude of 7^1=94. 5±0. 5 



km 



Therefore, the only parameter missing for an 
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accurate mass determination in X1822-371 is the radial 
velocity curve of the companion star. 

Hel absorption features attributed to the irradiated 
hemisphere of the companion star were presented in Har- 
laftis. Home & Charles (1997), Jonker et al. (2003) and 
Cowley et al. (2003). They measured K2 > 225 ± 23 
km s~\ > 280 ± 26 km s^^ and > 234 ± 20 km s'^ 
respectively by fitting sine waves to the Hel radial ve- 
locity curves. However, as Jonker et al. pointed out 
the spectroscopic Tq is delayed by O.lPorb with respect 
to the photometric eclipses and pulsar ephemeris and 
hence asymmetric irradiation of the Roche lobe needs to 
be invoked. 

In Casares et al. (2003; hereafter C03), the first un- 
ambiguous detection of the donor star is reported using 
the NIII A4640-Bowen emission line with a velocity am- 
plitude of Kem = 300 ± 8 km s^^. The phasing of the 
NIII velocity curve is perfectly consistent with the pulsar 
ephemeris with no need for ad hoc asymmetric illumina- 
tion effects. Furthermore, they show that the Hel absorp- 
tions come from the accretion stream and not from the 
companion star. This naturally explains the phase offset 
with the pulsar ephemeris observed in previous works. 

However, the NIII-Bowen emission line is excited 
on the inner hemisphere of the donor star and only 
provides a lower limit to the true ii'2-velocity. Narrow 
high-excitation emission lines originating from the com- 
panion star have also been detected in the prototypical 
LMXB Sco X-1 (Steeghs & Casares 2002), and in the 
black hole X-ray transient GX339-4 during its 2002 
outburst (Hynes et al. 2003). These narrow features 
are strongest in the Bowen blend, a set of resonant 
emission lines (mainly CHI and NIII A4640) excited by 
photoionization and the Bowen mechanism (McClintock, 
Canizares & Tarter 1975). The Roche lobe shaped donor 
star intercepts the energetic photons from the inner 
accretion disc resulting on the observed optical emission 
lines from its surface. The Doppler motion of these lines 
trace the orbit of the companion star and provide the 
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first reliable binary parameters in persistent LMXBs. 
However, the secondary is not uniformly illuminated by 
the ionizing photons and, therefore, the radial velocity 
curve of the sharp features traces the motion of the 
reprocessed light-center rather than the true center of 
mass of the star. 

In this paper we model the deviation between the re- 
processed light-center and the center of mass of a Roche 
lobe filling star (the so-called " K-correction" ) in a per- 
sistent LMXB, including screening effects by a flared ac- 
cretion disc. We apply our simulations to the case of 
X1822-371 and conclude that it must contain a massive 
NS with Mns ~ l.eiM© - 2.321/©, for any possible 
value of the disc flaring angle. 

2. THE "K-CORRECTION" FOR (EMISSION) LINES IN THE 
IRRADIATED FACE 

K-correction for radial velocity curves derived from ab- 
sorption lines formed on the companion's photosphere 
was presented by Wade & Horne (1988; hereafter WH88). 
In this case, absorption lines are quenched in the heated 
face and hence the observed K-velocity provides only an 
upper limit to the real K2, i.e. the true velocity of the 
center of mass of the companion star. Here, we study 
the K-correction for the case of irradiated-induced spec- 
tral lines originating from the heated face of the donor 
star, which is the opposite case to the one presented in 
WH88. This problem was flrstly addressed by Horne & 
Schneider(1989; hereafter HS89) but it was only applied 
to the particular case of nova V1500 Cyg without screen- 
ing effects by a possible accretion disc. Also Beuermann 
& Thomas (1990) studied this problem for the case of 
the cataclysmic variable IX Velorum. They considered 
only two different accretion disc angles together with 
the orbital parameters of this binary. In this paper we 
tackle the problem in a general approach and derive a 
K-correction expression for all possible q values and ob- 
scuring disc elevations. 

Following the prescription of HS89, the light-center of 
the emission line region(Zc) is given by Ic = a2 — d, where 
d represents the light-center displacement from the cen- 
ter of mass of the donor star and 02 = a/(l -|- g) is the 
distance between the center of mass of the companion 
and the center of mass of the system, q — M2/M1 is 
the mass ratio and a the binary separation. Using this 
nomenclature, the difference between the K2 and the ob- 
served emission line velocity {Kf,„i) is 

= K2-K,„, = K2^^ 

and therefore, 

K K 
l_rf(i±£) = l-/(l + g) 

a 

where / is an dimcnsionlcss factor which represents the 
displacement of the emission light-center from the center 
of mass of the companion in units of a. In the limit 
case where all the emission line is formed at the inner 
Lagrange point (Li) 

/ = ^ ^ 0.5 + 0.227 log g (Plavec & Kratochvfl 1964) 
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Fig. 1. — / factor or displacement of the emission light- 
center from the center of mass of the companion, as a function 
of q. Dashed lines represent the limit cases corresponding to 
/ ~ sin^ aM (lower curve) and / ~ (upper curve) . A more 
realistic upper limit to /, derived from the irradiation model is 
represented by a dotted line with the best polynomial fit overplot- 
ted in solid line as well. 

where 62 is the distance between Li and the center of 
mass of the companion star. On the other hand, we can 
consider the extreme case where the emission line comes 
from the irradiated point with maximum radial velocity. 
This point is located at the limb of the irradiated region 
and holds an opening angle um above the plane of the 
binary ^. This is the highest possible angle for the com- 
panion to be irradiated and it has been tabulated by R. 
Pennington in Table A3 of the Appendix to Pringle & 
Wade (1985). An approximate analytical expression, as 
a function of q, is provided by Paczynski equation (1971), 

sinaM = — = 0.462(^— )i/3. (2) 

a 1 + 9 

Then, if we assume that the Roche lobe is spherical 
(which is a good approximation when we are far from Li) 
it is easy to show that the minimum / displacement is 
given by ~ sin^ um ■ Therefore, in any possible situation 
/ will always be constrained between these two cases i.e. 

0.5 4- 0.227 log g > / > sin^ 

These limits are presented in Fig. ^ plotted in dashed 
line. Here we have used the exact aM values calculated 
by our code (see below) which obviously coincide with 
the ones tabulated by R. Pepringter. 

However, the upper limit obtained is purely theoret- 
ical since Li will never be the only irradiated point of 
the donor star. A more realistic upper limit can be ob- 
tained by simulating emission lines formed by isotropic 
irradiation on the inner hemisphere of the companion in 
a Roche lobe geometry. We have divided the Roche lobe 
flUing star in 10^ triangular tiles of equal surface. The 
flux from the X-ray source which is reprocessed in each 
tile will be given by, 

T-. Fx(l-e) cos Ox „ 
t^r ^ J3 1* 

where Fx is the X-ray flux emitted by the compact 
object, at the area of the tile, r the distance between the 

* Note that for this case to work the accretion disc must screen 
the rest of the inner face of the donor 
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Fig. 2. — Radial velocity curves for q = 0.6, i = 90° and four 
different disc obscuring angles a. The plot shows how the curves 
start to deviate from sinewaves for low a values. The amplitude 
of the velocity curves depends on a as well. Note that the velocity 
values have been rescaled by normalizing to v = 100 km s~^ the 
maximum of the radial velocity curve corresponding to o = 18° . 
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Fig. 3. — Radial velocity curves for q = 0.6, a = 0° and three 
different values of the inclination angle. The maxima of the radial 
velocity curves are closer to orbital phase 0.25 for the lower incli- 
nation cases. Note that the velocity values have been rescaled by 
normalizing to v = 100 km s~^ the maximum of the radial velocity 
curve corresponding to i = 90°. 



tile and the center of mass of the primary star, cos 9x the 
projection factor of the tile as seen by the X-ray source 
and e the albedo. The total reprocessed emission line 
flux seen by a distant observer is obtained by integrating 
Fr cos OyPg over the Roche lobe geometry, where cos^i, 
is the projection factor of each irradiated tile as seen by 
the observer and Pg represent a local gaussian profile. 
We have not included limb-darkening effects in the sim- 
ulation because these are expected to be negligible in 
the optically thin case that we are considering here. To 
test this hypothesis, we have also included in our sim- 
ulations a linearized limb darkening, appropriate to the 
continuum (i.e. optically thick case), and found it had 
no impact on the results. 

Radial velocity curves were subsequently computed by 
measuring the position of the maxima in the synthetic 
profiles (which are not symmetric), and if em is deter- 
mined from the maxima of the radial velocity curves. 
Therefore, for a given K2 eq (1) gives the / values as a 
function of q. These values, as expected, are not found to 
depend on either e nor Fx- The synthetic f{q) curve is 



plotted in Fig. ^in solid line and a 4th order polynomial 
fit (which is accurate to better than 1%) to this curve 
yields 

f{a = 0°) = 0.116 + 0.946q- 1.947(72 + 2.245q^ - 1.002g4 

Then, the real value of / will be constrained between 
these two limits i.e. /(a = 0°) > / > sin^ a^i and will 
depend on the size of the companion (i.e. the mass 
ratio) and the height of the accretion disc rim obscuring 
the heated star. 

In a second step, we have incorporated shadowing ef- 
fects by a flared and opaque accretion disc with an open- 
ing angle a independent of azimuth. Figure |3 presents 
the synthetic radial velocity curves computed for four 
different disc obscuring angles. The plot shows how the 
radial velocity curves start to deviate from pure sinewave 
functions when a decreases. Because the velocity am- 
plitude is strongly dependent on a so it does the K- 
correction. Therefore, using our binary code we have 
computed the relation between the ratio and q for 
different values of the flaring angle a and these are plot- 
ted in Fig. 0]in solid hue. 

In Fig. |5l we also explore the effect of the inclination 
angle in the shape of the radial velocity curve. This plot 
shows that the morphology of the radial velocity curves 
clearly depends on the inclination angle and hence, -^^^ 
could be a function of this parameter as well. Using our 
numerical model we have found that (in the most ex- 
treme cases when q > 0.6), -^^^ can be 10% smaller 
for low inclination angles than for i ~ 90°. To illustrate 
the effect of the inclination angle we have also computed 
the K-correction for i — 40° and different a values, and 
these are overplotted in Fig. ^in dash-dotted line. 
The synthetic K-corrections produced by the numeri- 
cal code are very smooth and can be well described by 
fourth-order polynomial fits. Table 1, lists the coeffi- 
cients of the fits for the two values of the inclination 
considered and different opening angles (from 0° to 18°) 
which cover the typical values observed in XRBs and cat- 
aclysmic variables. These fits are accurate to better than 
1% between q ~ 0.05 and q 0.83. For the case of in- 
termediate inclinations the interpolated solution between 
i = 90° and i = 40° is accurate to < 2%. For i < 40° the 
i = 40 value gives also a very good aproximation (< 2%). 
Note that the value of obtained for i = 90° yields 
a secure lower limit to the masses of the binary compo- 
nents for all possible cases. 

For intermediate a and i values one needs to interpolate 
the K-correction solutions and not the coefficients. 

Hence, for a given mass ratio q, the K-correction is 
provided by: 

^c^No + Niq + N^q^ + Nsq"" + N^q" (3) 

Note that, since the K-correction is only meaningful 
if a < au (otherwise the companion will be completely 
shadowed by the disc), eq. (2) sets a lower hmit to q for 
each a, (see Fig 4.). For this reason the first values of q 
considered in the fits for a= 10°, 12°, 14°, 16° and 18° 
are q— 0.1, 0.15, 0.26, 0.4 and 0.58 respectively. As a 
test, we have computed the K-correction for the binary 
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Fig. 4. — K-correction as function of q for different disc opening 
angles. Dotted lines represent the K-correction computed by our 
model for i = 90° and disc opening angles 0°, 8°, 14° and 18°. 
The solid lines mark the polynomial fits to the computed curves 
and the dash-dotted lines show the case of j = 40°. The dashed 
lines represents the same limits as in figure 1. 



Fig. 5. — Determination of the neutron star mass in X1822-371. 
The solid lines represent the results of the model for the open- 
ing angles considered (ol ~ 12° and a ~ 14°). The dashed lines 
represent the uncertainty in the orbital inclination. We have also 
marked in dot-dashed lines the measured K^m = 300 ± 15fcms~^. 
The possible neutron star mass range is restricted by the shadowed 
region. 



parameters of nova V1500 Cyg and obtain = 0.64 
in excellent agreement with HS89. 

3. THE "K-CORRECTION" FOR X1822-371 

Now we can apply our computed K-correction to the 
case of X1822-371 and derive the masses for the two com- 
ponents. In C03 a Kf.m= 300 ±8 km s~^ was reported 
through the detection of the NIII A4640-Bowen emission 
line arising from the companion star. This value, com- 
bined with the previous determination of the inclination 
angle and the pulsar's radial velocity curve, yields ab- 
solute lower hmits Mi > 1.14M0 and M2 > O.36M0, 
assuming K2 > Kem (see COB). 

On the other hand, the mass function equation can be 
expressed as 

Ml = —2 = 

2TrG sin"^ i 

where P is the orbital period and G is the gravitational 
constant. In the case of X1822-371, where i and Ki are 
known, this equation provides a direct relation between 
q and Mi. Therefore, this can be combined with our 
synthetic K-correction curves to yield Kem-Mi relations 
for different disc opening angles a. 

C03 obtained Kem — 300 ± 8 km s^^ by fitting the 
NIIIA4640 spot in the Doppler map with a 2D-gaussian 
function. Instead, we have computed the centroid of the 
spot intensity distribution and obtain a more conserva- 
tive value of Kem= 300 ±15 km s""'^, which we prefer 
since it does not impose any model to the fit. In order 
to test whether the K-correction derived by our model 
can be used to correct the Kem value extracted from the 
Doppler tomography we have simulated irradiated pro- 
files (using our derived parameters of X1822-371, see be- 
low) for the orbital phases of the C03 observations. These 
were smeared according to the exposure times used and 
subsequently degraded to the spectral resolution of the 
observations. Finally, we added noise to reproduce the 
observed spectra. A Doppler map was then computed 
and the centroid of the spot was found to be in excel- 
lent agreement, within the errors, with the value above. 
Furthermore, we note than the same K^m was obtained 
by SC02 in Sco X-1 from the radial velocity curve and 



Doppler tomography. 

Under the assumption that the Bowen emission lines are 
formed on the critical Roche lobe of the companion the 
detection of this emission at Kem — 300± 15 km s~^ di- 
rectly yields Mns > 1.35M0 and M2 > O.4M0 through 
the geometrical limit / > sin^ctM- However, realistic 
mass determinations require the knowledge of the disc 
rim elevation obscuring the companion. De Jong et al. 
(1996) derived a mean disc opening angle of a ~ 12° 
by studying the effects of X-ray reprocessing in LMXBs. 
This value was also derived by IIM89 for the highest disc 
structure (i.e. the bulge at phase '^0.8) in X1822-371 
by modeling simultaneous optical and X-ray EXOSAT 
lightcurves. However, in HNOl a higher angle of a ~ 14° 
is obtained by fitting X-ray RXTE-I-ASCA hghtcurves. 
These estimates correspond to the disc elevation at the 
phase of the bulge, and therefore, must be treated as 
upper limits to the height of the disc rim obscuring the 
companion, which is constrained to —0.04 < cj) < 0.04. 
In any case, these opening angles are expected to be, 
in general, only slightly higher than the real ones, and 
hence, we have calculated the i^em-Mi relation for both 
disc thicknesses as it is shown in Fig. O If we take into 
account, the opening angles considered, the inclination 
uncertainty and the errors in the observed Kem, Fig- ]M 
shows that the mass of the neutron star is constrained to 
the range 

I.6IM0 < Mats < 2.32M0. 

We would like to remark that this provides a conservative 
lower limit to Mjvs- If the accretion disc were thinner 
the mass would be even higher. For instance, if a = 0° 
(i.e. no disc), Mjvs must be ~ 3M0 in order to reproduce 
the observed Kem at ~ 300 km s~^. Besides, since the 
neutron star mass and the mass ratio are related by the 
mass function we obtain 0.24 <q< 0.27 which leads to 

O.44M0 < M2 < O.56M0. 

4. DISCUSSION 

The detection of irradiated-induced emission lines orig- 
inating from the heated face of the companion star can 
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yield important restrictions to the mass of the binary 
components. Kem contains intrinsic information about 
the companion's regions responsible for reprocessing the 
X-ray/UV radiation. Since the orbital parameters of 
X1822-371 are tightly constrained, the key missing in- 
formation for a realistic simulation are the structures 
present between the companion and the X-ray source. 
Our simple model only takes into account screening ef- 
fects by a flared accretion disc. However, the presence 
of a disc corona scattering the X-rays/EUV radiation 
might modify the results of our simulations since irradi- 
ating photons could reach (otherwise obscured) regions of 
the Roche lobe with lower velocities. Consequently, the 
scattered radiation is expected to push the K-correction 
to lower values and hence the NS to higher masses. As 
a test we have included emission from a 'corona' which 
has been roughly simulated by five X-ray spots symmet- 
rically distributed around the compact object at a dis- 
tance of 0.2a (i.e. the 'corona' radii derived in HNOl). 
The integrated flux of this simplifled 'corona' was scaled 
to be 1% of the accretion luminosity i.e. 10^^ erg s~^ for 
an assumed compact object luminosity of 10^* erg s— 1. 
Our simulation shows that the effect of such a 'corona' in 
the K-correction is completely negligible. Even if we rise 
the luminosity of the corona to 10% of the central X-ray 
source and take a > 10°, the effect in Kem is lower than 
1%. In this extreme case, the mass of the NS would only 
rise by ~ 0.03 M© with respect to neglecting the effect of 
corona scattering. In summary, our simulation provides 
compelling arguments for the presence of a massive NS 
in X1822-371 

To date, the most accurate neutron star mass deter- 
minations have been provided by studying radio pul- 
sars. These works reflect that the mass of the neu- 
tron star at formation is ~ 1.35 ± O.O4M0 (van Kerk- 
wijk 2001). However, recent studies have revealed ev- 
idences for heavier neutron stars in old white dwarf + 
pulsar systems, especially in J0751-I-1807 where a mass 
between 1.6-2.8 Mq has been reported by Nice, Splaver, 
& Stairs (2004). There is no compelling evidence for 
massive neutron stars in active XRBs yet, but possible 
candidates are Cyg X-2, with A4ms = 1-78 ± O.23M0 
(Casarcs ct al. 1998; Orosz & Kuulkers 1999) Vel X- 
1, with Mns = 1-86 ± 0.33Mq (Barziv et al. 2001) and 
V395 Car, with Mx = 1.9-4.1M0 (Shahbaz et al. 2004; 
Jonker et al. 2005). In the latter case the nature of 
the compact object is an open issue since NS signatures 
(X-ray burst, pulses) have not yet been detected. The 
evidence for a massive NS in X1822-371 reinforces the 
theory that LMXBs could harbour heavy neutron stars 
due to sustained episodes of high mass transfer. 

4.1. The nature of the companion star and evolutionary 
history of X1822-371 

Faulkner, Flannery & Warner (1972) showed that the 
mean density of a Roche lobe-flUing companion star is 
determined solely by the binary period P: 

p ^ 113P^-2 g cm-3 

In the case of X1822-371, with Phr = 5.57, we obtain 

p = 3.7 g cm""^. This value is consistent with a 0.5Mq 
MO main sequence(MS) star (Cox 2000) and it is in 



excellent agreement with our derived mass. We also 
note that a Q.5Mq MS star would not have evolved in 
a Hubble time. However, if we assume that the NS at 
formation was ~ 1.4M0 it must have accreted ~ O.5M0 
from the companion to reach its current ~ 1.9M0 and 
hence, it must have had at least > IMq at the onset of 
the mass transfer phase. 

Parmar et al. (2000; see also Baptista et al 2002) 
demonstrated that the orbital period of X1822-371 is 
increasing at P = 1.78 x 10~^°. As it was discussed in 
HNOl, this period variation can be used to determine 
the mass transfer rate from the companion M2, which, 
in a conservative case, equals the accretion rate onto the 
NS (Ml). The Ml value imphed by our ~ l.QM© NS is 
^ 3 — 6 X lO^^M0yr~^ ( comparable to the one obtained 
for an assumed 1.4M0 NS) and yields Lx ^ lO'^^erg 
s~^. Since X1822-371 is an ADC source we only see 
X-rays scattered by the "corona" and there are only 
rough estimates of the true luminosity. For instance, 
Jonker & Van der Khs (2001) derived 2 - 4 x 

]^q37 gj-g g-1 f]-Qj]2 the spin-up and constraints to the 
magnetic fleld. This, in turn, imply Mi ~ lO~^M0yr~^ 
i.e. one order of magnitude lower than predicted from 
P. However, this discrepancy could be reconciled if Mi 
were driven by secular changes and not reflect P (e.g. 
magnetic cycles in the companion ;Hellier et al. 1990). 
HNOl proposed that X1822-371 might have had a 
evolutionary history similar to Cyg X-2 (Podsialowski & 
Rappaport 2000 ; Kolb et al. 2000) where the system de- 
scend from an intermediate mass X-ray binary (IMXB) 
with g > 1. In this scenario the secondary must have lost 
more than ~ IMq from its outer layers in order to reach 
its current mass, and it should expose CNO processed 
material from the inner layers. Unfortunately, the 
companion star is not directly detected spectroscopically 
and hence the spectral type and the history evolution 
are unconstrained. However, the emission line spectra 
reported in C03 seem to show "normal-strength" CIII 
lines within the Bowen blend and also a strong OVI at 
A3811. Our observed masses, orbital period and Mi 
could be accommodated by some IMXB evolutionary 
sequences with an initial companion's mass and orbital 
period of ^ 3.5M0 and ^ 20hr respectively (Ph. 
Podsialowski private communication). The currently 
observed q ^ 0.25 implies that X1822-371 must have 
long passed the thermal mass transfer phase and its 
evolution should be currently driven by magnetic brak- 
ing phase with P < to sustain mass transfer, in clear 
contradiction with the observed P > 0. Alternatively, 
mass transfer could be sustained by X-ray irradiation 
as it was suggested in Pfahl, Rappaport & Podsialowski 
(2003). A fundamental test of the evolutionary history 
of X1822 will be to perform new high resolution optical 
and UV spectroscopy to search for evidence of CNO 
processed material (eg. XTE J1118-I-480, Haswell et al. 
2002). 

We thank Ph. Podsialowski for his useful comments. 
We are also grateful to the referee for his useful com- 
ments and suggestions which have improved an earlier 
version of the manuscript. JC acknowledges support by 
the Spanish MCYT grant AYA2002-0036. 
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TABLE 1 

K-OORRECTION FIT COEFFICIENTS 









= 90° 










= 40° 






a 


ATq 


Afi 


N2 


N3. 


N4, 






N2 


N3 




0° 


0.888 


-1.291 


1.541 


-1.895 


0.861 


0.886 


-1.132 


1.523 


-1.892 


0.867 


2° 


0.903 


-1.151 


1.591 


-1.914 


0.844 


0.897 


-1.159 


1.591 


-1.935 


0.860 


4° 


0.922 


-1.213 


1.855 


-2.358 


1.106 


0.904 


-1.477 


1.576 


-1.948 


0.888 


6° 


0.942 


-1.230 


1.877 


-2.336 


1.070 


0.929 


-1.313 


2.147 


-2.733 


1.264 


8° 


0.982 


-1.430 


2.552 


-3.284 


1.542 


0.975 


-1.559 


2.828 


-3.576 


1.651 


10° 


1.015 


-1.444 


2.426 


-2.957 


1.333 


1.015 


-1.590 


2.696 


-3.233 


1.439 


12° 


1.093 


-1.901 


3.593 


-4.280 


1.888 


1.130 


-2.160 


4.115 


-4.786 


2.062 


14° 


1.137 


-1.418 


1.324 


-0.679 


0.000 


1.165 


-1.617 


1.565 


-0.797 0.000 


16° 


1.326 


-1.871 


1.745 


-0.782 


0.000 


1.428 


-2.362 


2.421 


-1.122 


0.000 


18° 


1.496 


-1.645 


0.665 


0.000 


0.000 


1.579 


-1.859 


0.767 


0.000 


0.000 
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